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Oa.tbA  Qieoxy  of  Nonst&tionaxsr  ScatterASN 


I.'  A.  KLOOSBIN 

Ukr*  Zis*  Zh.^  v6I.  5,  No.  $,  i960,  pp.  6S0  -  627 


A  strong  shockvave-f root  vbicli  moves  in  the  atmospbere  of  tho  earth  ( or  even  la 
a  starxjr  savsloips  or  Is  interstellar  sipace)  raAtates  light;  the  radiant  energy  Is 
here  ahsoxhed  hy  the  gas  ahead  of  the  vavefrent  to  some  degree.  Radiant  energy 
frcoi  a  strong  tfhookwsve  fkont  considerably  affects  the  magnitude  of  the  ten^rature 
at  the  mnrefront  and  detenslnes  Its  eonflguratlon.  ^  The  effect  of  Inter¬ 
ference  from  radiation  oa  the  eharaoterlsties  of  strong  "terrestrial”  shockwaves 
vas  eonsidered  by  ZA.  fi.  Zel'dovlch  [IJ  and  by  XU.  R.  Raiser  [2].  The  effect  of 
radiation  ftom  a  dhodcwave  front  vhleh  moves  In  the  shell  of  a  star,  on  the  magni¬ 
tude  of  the  tenqperature  Juiqp  at  the  wavefront  Is  analyzed  ih  [3].  Uhe  results  of 
the  abovo  vork  oonfixm  the  aeeesslty  of  taking  into  account~i^he  Interaction  of 
strong  shockwaves  with  radiation  indejiendentjiy  of  whether  they  move  In  the  atmos- 
Jiiere  of  the  earth,  a  star  or  In  Interstellar  space. 

But  up  to  reoently^  it  vas  eustooary  to  solve  a  system  of  equations  of  motion 
together  vltdi  the  equations  of  radiation  transfer  la  Investigations  of  these  pro- 
'  blems.  On  the  vhols«  this  ^vtem  is  coopUcated  and  Is  solved  only  by  numerical 
methods.  Undoubtedly  the  use  of  ccmteoporaxy  methods  of  the  theory  of  light  scat¬ 
tering  eaa  sffoxd  new  possibilities  in  this  direction.  Put  only  the  case  of  light 
scattering  in  a  medium  whose  boundaries  remain  fixed  has  been  solved  In  problems 
of  this  theoxy.  In  the  case  of  light  scattering  ahead  of  or  behind  a  shockwave 
front,  tho  boundary  of  the  medium  which  Is  scattered  moves  relative  to  the  gas. 

Bsnoe,  before  tgplylng  the  modern  theory  of  light  scattering  to  g^sdynaalcs  probleau, 
it  Is  necessary  Uf  develqp  a  theory  of  nonstationary  light  scattering  in  a  medium 
with  moving  boundaries  • 
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VBmSSSAS*  8SATEMEHIS  OF  THE  THEORY  OF  ROHBTATIOIIARY  LIGET  SCATTBRIHQ 
Of  all  tbe  Mtho&s  Of  sestteriaK  theory^  the  nsthod  of  detexninlng  the 

;ipro*batKULll^  that  a  Qjaantan  votUfl  leave  a  jnedlim  vhleh  vaa  projposed  V.  V.  Sobolev 
[4J  eeeoed  ttoat  anltabla  to  ue*  IVm  this  method,  the  probability  p(t^T)  Is 
found  that  a  Qiastm  ahsoorbed  to  the  ogptlcal  depth  t  at  zero  time  vlU 

leave  vlthln  the  time  t  •  Two  cases  hence  exist:  l)  Most  of  the  tljas  the 
Is  ia  the  ahsortMd  statej  2)  most  of  the  time  the  qiumtimi  resides  on  the  path  het» 
veea  too  ahsotptions*  Sie  first  easa  has  meaning  la  the  atmosphere  of  the  earth 
a&d  In  tile  onvelcpes  of  stars^  tile  eeeoad  In  interstellar  space  and  In  the  tqper 
layers  of  the  atmoaiihiere  of  stars* 

7*  V.  Sobolev  formulated  Integral  and  equivalent  differential  equations  for  the 

probahUtty  that  a  quanttm  vould  ^clt  Aran  the  medium  In  both  cases*  Xn  particular^ 

the  aolutiou  for  the  function  pCf^u)  ,  vhere  u  «>  is  nondlmenslonal  tlmft>  Is 

^1 

•  -(1 - 2^)u 

U)  1  r (x  cos  xt*alaXT)e  ^  ^  -g  , 

to  the  first  case^  >diere  X  le  the  probtiilllty  of  the  etrrvival  of  the  quanbun 
after  the  act  of  ahaoiptloni  t^  the  mean  time  spent  by  the  qjuanttmi  on  the  path 
hetveen  tuo  aeatterlnfi  acta* 

Tor  pcaotlcal  ptorposea  the  determlBation  of  the  funetion  p(t)  «  the  probahll* 
itgr  that  a  Ugtt  quantun  absoxhed  at  the  optical  depth  x  vill  generally  leave 
the  meditmif  and  of  the  function  E(t)  ,  vhlch  detezmlaea  the  mean  time  the  quantum 
apands  la  the  sKdltaa,  is  of  oonslderohle  interest*  3hese  fonstlons  are  determined 
from  tha  expressions 


w 

PCt)  t(XtV)dVit 
0 
«0 

Z(t)  m  J  up(T/n)dn. 


V*  Sdbelev  found  Sx^gsaX  and  equival^xb  differential  equations  for  PCv) 
sad  I^t)  aa  well  as  the  solutions  for  these  functions  in  the  case  of  a  bounded 
sod  for  X  ■  1  .  Bte  function  P(t)  for  a  seml^infinite  medtum  is  fouad 


from 


(4)  pCt}  -  (1  -  -  X)e" 

Let  us  note  that  for  the  case  of  a  seal-infinite  nedlim  ve  have  found  the  solu¬ 
tion  of  the  apiprpprlate  dlffsreatl-al  eQ]uatlen  for  Z(t)  as 

(5)  ^  «  (1  -  X)Z  -  XP 

vlth  tile  houndazy  eondltloas  Z^q^vZXO)-*-  X* 
fhla  soJjtttloa  can  be  tfrltten  thus 

(6)  Z(t)-— ~=r-  U  ♦  (1-Vl  -  X)Tle 

2Vi^ 

Xn  the  ei^resslons  presestedf  x  •  knx  ie  the  optical  depths  there  k  is  tiie 
absorption  coefficient  per  particle j  a  the  nunflier  of  particles  per  unit  voLoae 
and  X  a  gecmetrle  coordinate. 

Shis  veil-behaved  and  coopletely  eonelvsive  theory  can  also  be  called  to  the 

case  of  light  scattering  vith  a  movlsg  bottndaxy,  Xo  this  case^  the  qptlcal  depth 

will  he  considered  as  the  position  of  the  honndax;^  at  a  given  Instant*  Let  us 

laaglne  the  medium  to  be  semi- infinite  (in  Optical  thlelsaess),  ItC*^  bounded  on 

4t 

one  side  hX  ^  boxuidary  which  moves  with  the.  constant  velocity  ^ 

or  out  of  the  medlimi.  Hence,  let  tts  imagine  that  the  light  la  not  scattered  on 
'the  other  side  of  the  moving  boundary*  Por  simplicity  of  the  computatlonSj  let 
US  be  limited  to  the  case  of  a  one -dimensional  medium.  Let  US  also  Imaglpe  that 
the  probabUlV  of  survival  of  a  quantum  after  scattering  X  la  Independent  of 
tbe  optical  thickness  and  that  the  probability  of  scattering  on  both  sides  is 
the  same* 

XROBABIim  TEAT  A  QUANTUM  WILL  LEAVE  A  MEDIUM  WPES  A  HOVIIK}  bOUHDARI 

AT  THE  SIDE  OP  THE  MEDIUM 

A  boundary  movliig  with  velocity  v  will  traverse  a  path  uv  vtthln  the  tima 
tt  *  Seace^  In  this  case  it  is  required  to  search  separately  for  tho  probabilities 
that  a  qjoantum  will  leave  the  optical  depths  r  <  uv  and  T>uv'  *  If  r  b-  xiv  , 
then  the  probability  that  the  quantum  will  leave  vithottt  scattering  vlthln  the 

time  u  vlU  equal  ^"(t“w)-u  probability  that  the  quantum  will  leave 

• 

after  scattering  can  be  found  by  multiplying  the  probability  of  the  transition  of 
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a  quantum  from  the  depth  t  to  the  depth  t'  vlthln  the  time  u'  ,  vhlch  Is 

|-u*  probability  that  the  quantum  will  leave  the  depth  t’  at 

the  time  u  -  u*  -  pCt'*u*v^u>u' }  •  It  Is  hence  required  to  take  Into  considera¬ 
tion  that  actually  the  distance  fxca.  the. quantum  to  the  moving  boundazy  at  time 
«*  viUL  be  T*  -  u’v  .  Integration  should  be  over  all  t  froo  u'v  to  infinity. 

Bence^  making  the  change  of  variable  t**  »  t'  -  u'v  ,  we  find  that  for  all  T>ttr 

U  ^ 

(7)  .  pCt.n)-|  1/.-“’  du*  f 

O  o 

The  ejgparessloa  for  pC^iu)  for  the  optical  depth  T<ttr  has  a  less  eomplex 
form.  She  probability  that  a  quemtum  vlll  exit  without  scattering  Is  ~|e~'^  •  The 
probability  that  a  qoantimi  will  exit  after  scattering  la  composed  of  two  factors: 
the  first  detezmlnes  the  parobahlll^  that  a  quant \mi  will  exit  after  ssatterlng 
with  transitions  from  the  optical  depth  X  to  the  d^pth  t*  dinrlng  the  time  the 
movable  boimdary  has  not  yet  achieved. the  eptleal  depth  x  f  i.e.,  In  the  time 

the  second  factor  takes  into  moeount  the  probabillly  that  a  U^t  quantum 
will  exit  from  the  meditta  which  has  remaiaisd  In  the  absorbed  state  a  time  greater 
than  ^  t  has  been  carjpled  out  of  the  medluia  along  with  an  atom,  then  reradlated  to 
the  side  of  the  aedium  and  has  been  aheosbed  at  the  optical  depth  t*  .  Ihus,  the 
iategsal  eqcatloa  for  pCt^u)  for  all  t  <w  is 


(8) 


u  « 

J' a'*  p(T*,u-u'idT' . 


jr  o 

V 

She  appropriate  differential  equation  and  boundazy  conditions  can  be  obtained 
froa  (7)  and  (8).  The  dlfferentleuL  eqic&tlon  can  be  solved  by  sepcuratlon  of 
taxiablcs.  Bowever^  It  should  be  noted  that  these  equations  are  vezy  awkward  and 
the  solutions  obtolned  have  a  feem  Buch  that  it  Would  be  difficult  to  obtain  d 
Clear  pbyaieaX  picture  and.  estimate  of  the  quantities  for  comparison  with  the 
results  of  observations.  Hence,  instead  of  the  solved  equations  (7)  and  (8),  it 


Is  eaqpedlent  to  obtain  appropriate  equations  from  then  to  deteznlne  the  prohahlll'ty 
F(t)  that  a  q.\ian;i:\Bi  vould  leave  the  aedlm  at  ecias  Instant  and  Z(f)  the  average 
tine  the  quantum  stays  in  the  medium  after  ist  to  the  depth  t  ,  vhlch  are  defined 
hy  (2)  and  (3). 

In  the  case  ve  are  considering  t 

is)  p(f,tt)du  -yp^ft#n)dtt  *  J pg(t,u)dt^, 

o  o  ^ 

vhere  Is  the  e^qpresslen  C7)>  tdiich  defines  the  porfibabllltgr  cooBldoMd  at 

tt<^  and  for  u>^  Substituting  (7)  and  (8)  late  (9)s  ve  oibtaiA  after 

Integration  vlth  recgpect  to  u 


PC^Odt’  * 


Prom  this  equation^  ve  have  ^  dliXareatlal  equation  determining  P(t)  t 

Besce^  It  folXovs  that  the  funotioa  ?(t)  ecu  he  vrltten  as 

E-fcf 
“i*  * ' 

vhere  Ibhe  are  eonstants  detexmlaed  trcn  tSw  hoxindaxy  copSitlGocif  the  roots 
of  the  characterlstle  eqpiatiea 

(13)  Ic* .  . 

Of  the  three  roots  of  this  eq;aati(»i,  one  Is  negative^  it  la  discarded  since  the 
ablution  must  remain  hounded  as  x  -*m  . 

Substituting  (12)  Into  (lO),  the  Integration  can  be  perfomBd,  the 
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coefficients  of  the  various  esg^onents  eqjuated  to  zero  and  conditions  obtained  to 
detezwLne  the  constants  c^  : 

i  ^  i  ^  i  i  V 

FOzBUlas  (12)  -  (lU)  cosQtletely  detemine  the  solution  of  our  problem;  th\is 

(«)  p(,).  -  V>.  .-v  ^  /v 

(JCj-l^Xl+k^k^v)  (kj^-k2)(l +k^k2v)‘ 

vhere  the  found  from  the  characterlstle  eq^uatien  equal: 


2  Jl+3v  cos 


/iS  +  i 

13  3 


An  aigproxliaate  fomula  for  P(r)  can  be  ebtedned.  Thiis,  for  v  4tl 

-  (1.X)V|  k^m  +  ; 

)  vf  -  -  / - 

p(t1  r  .  ve  +  (1 -ylTx)e”>^^  . 


and  for  v  ^  1 


1  +  Vl-A 


k,  •  1  +  ^;  k_  ; 

1  2v'  ^  V  '  . 


rH' 

SiBlXarly,  vs  find  the  Integral  equation  for  Z(t) 


«  V 

(19)  Z(t)  ■  y'\ip(T,u)du  - y'pj^(T,vOu  du+y^ P2(t,u)u  du. 

O  O  T 

Substituting  the  approadLaate  values  Pj^(t,u)  axid  ^  P2('<^>u)  according  to  (?) 
and  (8)  and  integrating  vltb  respect  to  u  ,  ve  find 


K-2^ 


7 


2(1- v)' 


z(t)  .  — I  •  ^  + 

2(1- t)2  *  (IWT)® 


f  e'(''‘’^P(t' )dt*  + - i — 5  /e'^'''‘''^P(T‘)4t'  - 

2(1+ vr  ir 

T 


(20) 


8^v 


T  T^T 
V 


rj^T 

V 


u-^f 


Je  ^  PCtOdT* - ^  J'  (-r-T»)e  ^  P(t' )dT'  + 


7  (t-tO 


5J5^  /.-<»*-'5z(t.)1t'  -  ^/e  ■  Z(,.)4t'  ♦ 


T  ~  ■  O 

•  i  *"  .x* 
o 

Xt  is  csfi/  to  oT^taln  the  differential  equation  for  £(r)  frm  this  equation 

frcm  idildh  it  Is  seen  that  the  solotlen  Xor  tbs  Ametlon  Z(t)  osa  he  rtpreseateft 
In  the  form  .  *k  T 

J22J  Kt) ♦  VJe 

Substituting  (22)  and  (12)  into  (20)  sai  integrstlngf  vs  find 


(23) 


B, 


3k^  *  2kj^  -  V 

2v3(k^  -  if  -  (1  *  v)(k^v  -  if 


B. 


^  (kj^  -  l)(k^v  -  iXd  +  ^ 

13ie  solution  for  Z(t)  Is  now  vrltten  as 
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Z(t)  ■  ^ 


kj,(k^+l)(lc^-l)(2v3(k^-l)2  .  (i+v)(k^y  -  1)^1 


(24) 


X  1+ 


(k^-l)(kj^v-l)[(l4>v)(k^v-l)-2v"^(k^-l)] 


2v3(k^-l)2  -  (I4v)(k^v-1)2 

viiere  m-S  for  1»1  and  m-l  for  i«  3  j  aiii  kg  are  determined  by  (l6). 
Expression  (24)  Is  considerably  slnpllfled  for  v>l  cmd  reduces  to 

1-X-r 


(25) 


PROBABILITY  THAT  A  QUAHTUM  Will.  LEAVE  THE  MEDIUM  FOR  MOTION  OP 

THE  MEDIIM  BOUNDARY 

Stairtlng  from  considerations  aimlogous  to  those  presented  above,  we  find  for 
the  probability  that  a  quantum  will  leave  the  medlxim  if  the  boundary  of  the  medium 
moves 

u 


(26) 


p(t,u)  .  I  *  I  y'e-“'  du'  y'e-l^-’’''^'’'lp(T",u-u-)dT”, 


o  0 

ai^  Integrating  this  expression  with  respect  to  u  ,  we  obtedn  an  Integred.  ejq>re8- 
sion  for  the  probability  P(t) 

T 


(27) 


P(t' )dT' . 


*  /.-<^'-^>P(T-)dr.  - 

The  differential  equation  for  P(t)  in  the  case  idien  the  boundary  of  the  medlTm 
moves  is 

■  (28) 


d^P  1  d^P  ^  ^  Ij^  Q 

V 


dT3  ^dT^ 

This  equation  is  easily  obtained  from  (27)  ai]d  also  from  (U)  by  changing  the  sign 
of  the  velocity  v  .  Representing  the  function  P(t)  in  the  form  of  (12),  we  find 
the  characteristic  equation  ffom  (28)  .  We  find  from  the  condition  of  bovindedness 
of  the  solution  as  t  -»  «  that  0^*0,  0,  so  that  k^<0  and  k^^O.  Fur¬ 

thermore,  auibstitutlng  (12)  into  (27)>  ve  find  the  condition  determining  the 


eoiurtant  Cg : 


(29) 


(1  - 


Hsace,  the  expreeslon  for  P(t)  In  this  case  is 

(30) 

vbere 


-kgT 

P(t)  -  (1-  kg)®  / 


(31) 


^  [a  y/T^co. -  I  =o.  i  -  i]  . 


(32) 


The  asproxtmate  fozmilas  for  v^l  are 

kj  - /iTj;  +  < 


P(t)  »  (1  -  Vl  ^)e 


(33)  • 

Similarly,  ve  also  find  an  Integral  equation  for  Z(t).  Substituting  (26)  into 

(3)i  we  obtain 


(3h) 


Z(t) 


2(l+v)  2(l+v) 


y^“(T-T*  )p(^,  + 


2(l-v)‘ 


r e-('‘'-'')r(,-)dT'  -  g/ 
_  (t»-t) 


»  (t'-t) 


P(t' )dT' 


-  y(T'-T)e  ^  P(T')dT«  +  5^^y"e‘^'^"‘''^Z(T')dT'  + 

1-  V  Y  T 


T  T 

F^gb  this  expression,  ve  obtain  the  differential  equation 


(35) 


d^Z 

dT3 


1  d^Z  §Z.  1^  7.  1  fp  ^  ^  1 
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Sabstltutlog  the  solution  for  Z(t) 
flni  that  the  coefficients  and 


here  in  the  form  of  (22),  as  before,  ve 
satisfy  the  conditions 


B. 


(36) 


'1  “  2  i 

^  3ki'<r+2k^-v  ^ 

„  1  B 

1  1  -  ®1 


Taking  into  account  that  c^  «  c^  s  o  in  the  case  idien  the  boundary  moves  out¬ 
wards,  ve  obtain  the  solution  for  Z(t)  as 


(37) 


Z(t). 


[l+(l-k2)T]e 


3k^v+2k2-v 


which  8l]iQ>llfles  considerably  for  v:^l  and  can  be  written  as 

(38)  Z(t)  >-^[1  +  2^  t]  e-\ 

It  is  easy  to  see  that  each  of  the  expressions  and  solutions  obtained  above 
transform  for  v  =«  0  Into  the  known  solutions  found  by  V.  V.  Sobolev. 

The  question  of  applying  the  theory  presented  to  practical  problems  related  to 
the  expansion  of  shockwaves  In  the  terrestrial  atmosphere  will  be  examined 
separately. 

Let  xis  note  that  the  equations  and  solutions  presented  here  were  partially 
published  In  [5]. 

In  conclusion,  the  author  Is  grateful  to  S.  A.  Kaplan  for  guidance  during  the 
research. 

L‘vov  Iftxiv.  March  I6,  i960 
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